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I. INTRODUCTION 


The investigation of the behavior of nuclei under extreme conditions has become an important tool to reveal new facets 
of nuclear matter. In particular, nuclei at the frontiers of the valley of stability constitute nowadays the most active 
research area of nuclear structure physics. With the advent of new radioactive beam facilities it is now possible to 
produce exotic nuclei that may have occurred naturally in the interior of exploding supernovas jjj. In short, extremely 
proton- and neutron-rich nuclei are now within experimental reach. Of specific interest to the present paper are the 
considerable experimental efforts to study nuclei with roughly the same number of neutrons and protons, N ss Z. 

On the theoretical side, the challenge is to investigate whether models, developed for—and using the phenomenology 
of—stable nuclei, can still be applied in these new, as yet uncharted regions and, if not, to propose new approaches 
to do so based on the data available up to now. One of the main open questions is the validity of the nuclear shell 
model with its traditional magic numbers and of the usual treatment of the residual interaction among the valence 
nucleons ]|||. 

The nuclear mass is a property of quintessential importance as it directly determines the stability of a nucleus. 
There are several theoretical approaches that reproduce the systematics of masses of nuclei and it is worthwhile 
to mention here two of them. The Extended Thomas-Fermi plus Strutinsky Integral [Q (ETFSI) is a high-speed 
approximation to the Hartree-Fock (HF) method with pairing correlations taken into account through BCS theory. 
In earlier versions a Wigner term was not included and this has been claimed to be the reason for the systematic 
calculated underbinding by about 2 MeV for even-even N = Z nuclei [||. This effect persists for iV = Z odd-odd 
systems and for N = Z ± 1 odd-mass nuclei but with less prominence. The mass formula based on the Finite Range 
Droplet Model (FRDM) || starts from a sophisticated liquid drop mass formula to which microscopic corrections due 
to shell effects are added. Both approaches have comparable numbers of parameters (about 15) and make reliable 
predictions with impressive success. In the FRDM and also in a recent ETFSI calculation jfj a Wigner (correction) 
term is included that specifically deals with the peculiar behavior of binding energies of N ss Z nuclei and has a 
cusp-like behavior for N = Z. This treatment is effective for known masses but, as the correction is entirely ad hoc, 
it has the drawback that an extrapolation to unknown nuclei can be dangerous. It is therefore of interest to develop 
models based on simple physical principles that can account for the behavior of nuclear masses at the IV « Z line. 

Many models have been used over the past years to investigate the structure of heavier IV s=s Z. We mention in 
particular recent applications of the Hartee-Fock-Bogolyubov (HFB) method that includes proton-neutron pairing 
correlations ||. This approach is tailor-made for the treatment of N « Z nuclei but has the drawback of the 
lack of particle-number projection. Shell-model calculations are generally extremely successful in reproducing 
spectroscopic nuclear data but require large configuration space diagonalizations. In addition, no reliable theoretical 
procedure exists for deriving the monopole part of the residual interaction which hinders a comparison with measured 
nuclear binding energies. An algebraic approach [Q, which has affinities with the one presented here, utilizes the 
concept of dynamical supersymmetry for the calculation of the binding energies in the sd shell but does not go beyond 
it. 

In this paper, the Interacting Boson Model (IBM) |ll|] in its isospin invariant version is applied to proton-rich 
N « Z nuclei. The goal is to obtain reliable estimates of binding energies of T = 0 and T = 1 ground states in 
self-conjugate (IV = Z) nuclei based on the concept of dynamical symmetry. The mass formula proposed is relatively 
simple and contains terms with an intuitively understandable significance. A particular ingredient is its treatment of 
the competition between isoscalar and isovector pairing. 


II. AN IBM-4 MASS FORMULA 

The Interacting Boson Model in its original version (IBM-1) [|l2| successfully describes collective aspects of nuclei 
through the use of s and d bosons which are thought to approximate pairs of valence nucleons coupled to angular 
momentum 0 and 2. No distinction is made between neutron and proton bosons. Whenever the difference between 
the neutron and proton fluids is thought to play a role, one is forced to use more elaborate versions of the IBM. 
The neutron-proton interacting boson model, or IBM-2, was introduced mainly to provide a microscopic foundation 
to the model [ p~3| ]. It uses as building blocks s and d bosons constructed from neutron-neutron (nn) and proton- 
proton (pp) pairs solely. In the third and fourth versions of IBM, IBM-3 and IBM-4, the isospin quantum number is 
introduced in a natural way. In IBM-3 the entire isospin triplet T = 1 is included, leading to nn, np, and pp pairs 
with T z = Tl, 0, —1 |l4j. The IBM-4 considers both T = 0 and T = 1 pairs; the T = 1 bosons are assigned an 
intrinsic spin S = 0 while T = 0 bosons carry an intrinsic spin S = 1 |li|. A justification of this choice is that the 
two-particle isospin-spin combinations ( TS) = (10) and ( TS ) = (01) are lowest in energy and that they give rise to 
an SU(4) algebra which is the boson equivalent of Wigner’s supermultiplet algebra [jl6| . 
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The mass region 28 < N w Z < 50 has a very rich structural behavior, presenting many aspects of nuclear motion. 
It is an ideal testing ground for various models since a proper description of the data relies on the interplay between 
T = 0 as well as T = I pairing and deformation-driving interactions. In addition, it is a region of intense experimental 
studies but with few experimental data available up to now. 

Very recently, the IBM-4 was applied to the spectroscopy of exotic N « Z nuclei in the pf' 5 / 299 /2 shell |Q. In this 
approach the IBM-4 Hamiltonian is derived from a realistic shell-model Hamiltonian through a mapping carried out 
for A = 58 and 60 nuclei. The boson energies and the boson-boson interactions are thus derived microscopically and 
no parameter enters the calculation (since the shell-model interaction is considered as an input). This microscopically 
derived Hamiltonian gives good results in 62 Ga (when compared to the shell model) and predicts the energy spectra 
of heavier N = Z nuclei (such 66 As and 70 Br). The approach is reasonably successful in obtaining a spectroscopy of 
low-spin states in N « Z nuclei. It makes use, however, of a complicated Hamiltonian and, moreover, calculations 
beyond 70 Br seem difficult. 

Prompted by these considerations, in particular the need for reliable binding energy predictions at the N = Z line 
and the existence of a microscopically derived IBM-4 Hamiltonian, we propose here a simple mass formula in the 
context of IBM-4. In previous work [Q one of us introduced an algebraic Hamiltonian (which can be regarded as 
the s-boson channel of the general IBM-4 Hamiltonian of Ref. |Q]) with the specific aim to study the competition 
between the isovector and isoscalar pairing modes in self-conjugate nuclei. The s bosons do give rise only to the 
symmetric representation of Ul( 6) but this approximation can be justified in even-even and odd-odd N = Z nuclei 
(the only ones considered here) where this is indeed the favored U(6) representation [^ 5 ). 

These previous studies |lh,h8[ suggest that the relevant terms in a simple IBM-4 mass formula must be taken from 
two different symmetry classifications: 

U(36) D (U l (6) D •••=> SO l (3)) 

® (uts( 6) D { ® u](3) } D S0 ^< 3 > ® S °s( 3 >) ■ m 

Due to the overall symmetry of the boson state under U(36), the Ul( 6) and Uts( 6) representations have identical 
labels [Ni ,..., N 6 ]. 

The terms chosen from (]l|) for the calculation of binding energies of N = Z nuclei are the following. First, the 
linear and quadratic Casimir operators of Us (6) are included. The total number of bosons, N, labels the symmetric 
representations of Ul( 6) and, as a result, these Casimir operators take account of the bulk properties of the nucleus 
and lead to a smooth variation of the mass with particle number. The next two terms to be included are the quadratic 
Casimir operator of SUts( 4) and the linear Casimir operator of Us(3). In the L = 0 channel considered here, these 
include the effects of T = 0 and T = 1 pairing interaction; SUts( 4) implies equal T = 0 and T = 1 strengths while 
Ug(3) breaks this degeneracy by removing all non-scalar Us(3) states from the low-energy spectrum. In Ref. JlJ] 
the quadratic Casimir operator of SUs(3) was considered; here the linear Casimir operator of Us(3) is preferred 
because it is shown since to be connected with the spin-orbit term in the mean field which breaks the isoscalar- 
isovector degeneracy in favor of isovector states |juj|| . The final term considered is the quadratic Casimir operator of 
SOt(3), T(T+ 1), which is known to represent the nuclear symmetry and Wigner energies. In summary, the following 
Hamiltonian is taken: 

H = E 0 + aC 1 [U L (6)] +/3C 2 [U L (6)] + 7 4[SU TS (4)] +£Ci[U s (3)] +rA[SO r (3)], (2) 

where Eq is the energy of the doubly magic core, specific for a given mass region. Note the omission from (|2|) of 
operators associated with Urs(6), Ur(3), and SOs(3); these are not needed because their effect is equivalent to the 
corresponding operators of Ul( 6), Us(3), and SOt(3). 

All operators in (|2]) mutually commute, except for C2 [SUts( 4)] and C\ [Us(3)] and hence the solution of H involves 
a numerical diagonalisation which is most conveniently done in the second basis in (|l|), labeled as [N] (\t)T <g> (As)S'). 
The matrix elements of C^SUtsG^)] can be calculated analytically |il|, 

V x N T l S s y T y s = ([N](\ t )T® (As)S|4[SUsr(4)]|[iV](A^)T® (X' s )S), 

V x N t T x S sXt x s = 2X tXs + 3 N + T{T + 1) + S(S + 1), 

V x n t t x S sXt . 2Xs+2 = [(A t - T)(A t + T + 1)(A S - 5 + 2)(A S + S + 3)] 1 / 2 , 

V x n t t x S sXt+2Xs -2 = [(At - T + 2)(At + T + 3)(A S - 5)(A S + S + l)] 1 / 2 , (3) 

while the other operators are diagonal with eigenvalues given by 
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([N}(\ t )T®(\s)S\C 1 [\J l (6)}\[N](\ t )T®(\s)S) = TV, 
([TV](A t )T ® (X s )S\C 2 \U l ( 6)]\[N](X t )T ® (A s )S) = TV(TV + 5), 
([TV] (A t )T <g> (X s )S\d 1 [U s (3)] | [TV] (A t )T ® (A s )5) = A s , 
([TV](A t )T ® (A s )5|C 2 [SOt(3)]|[A^](At)T ® (A s )5) = T(T + 1). 
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III. RESULTS 


We first apply the mass formula to TV = Z nuclei in the sd shell, from 16 0 to 40 Ca, where the experimental masses are 
well known Q]. The five parameters of the Hamiltonian (||) are adjusted to the binding energies of the lowest states 
with T = 0 and those with T = 1 of all even-even and odd-odd self-conjugate sd-shell nuclei. There are thus two 
data points per nucleus, which is crucial for a reliable determination of the parameters 7 , £, and rj. Binding energies 
are corrected for Coulomb effects according to the prescription given in j(|. A drawback of the present formula is the 
occurrence of a discontinuity at mid-shell which is related to a change of core (from 16 O to 40 Ca). To avoid these 
mid-shells effects, two different fits are performed for each half of the shell, a first one for nuclei from 18 F up to 28 Si 
(TV = 6 bosons) and a second one (with 40 Ca as a core) for nuclei from 38 K down to 30 P. In Table | the two parameter 
sets are given in the lines labeled il6 0 to 28 Si’ and i30 P to 40 Ca’. The major difference between the two sets is the sign 
change in a which is required since in the first half TV counts the pairs of nucleons added to 16 O while in the second 
half it counts the pairs subtracted from 40 Ca. One also notes that a has a larger absolute value in first half than in 
the second: this must be so since, in the sd shell, the binding energy per nucleon increases as the size of the nucleus 
grows. Furthermore, the parameters 7 , £, and ?? decrease (in absolute value), as a result of the average interaction 
strength which decreases with mass. Nevertheless, one notes that this decrease is stronger for 7 than it is for £, that 
is, the ratio |£/ 7 | is larger in the second half of the sd shell than it is in the first. Again, this is understandable 
intuitively because one expects the Wigner SU(4) symmetry to be increasingly broken by the spin-orbit term [Us (3)]. 
The resulting binding energies for each half of the shell are shown in Tables [0 and III, respectively. Also the isospin of 
each state is indicated as well as the difference A E between the calculated and measured binding energies. In Fig. |l] 
the differences in energy between the T = 1 and T = 0 states are compared with the observed ones and also with the 
semi-empirical formula for this quantity given in Ref. |f22[ . The root-mean-squared (rms) deviation is 0.876 MeV in 
the first half of the sd shell and 0.245 MeV in the second half. Since reasonable results are obtained with parameters 
that can be qualitatively understood from simple arguments, an extension towards the 28-50 shell can be considered. 

We begin with a discussion of the first half of the 28-50 shell, for nuclei ranging from 58 Cu to ' 8 Y. The ground 
state of all these self-conjugate nuclei is a 0 + , with either T = 0 in even-even or T = 1 in odd-odd nuclei, with the 
exception of 58 Cu which has a (J 7I ’,T) = (l + ,0) ground state. Up to 64 Ge the masses are well known and can be 
taken from the compilation of Audi and Wapstra jnj. Of the heavier TV = Z nuclei, also the masses of 72 Kr and 
74 Rb are listed by Audi and Wapstra. The masses of° 6 As and 68 Se are available from a recent measurement Q and 
that of 76 Sr from Ref. j24j[. The latter experiment also gives a mass for 68 Se but since it is far off the systematics 
of Audi and Wapstra, the result from [^ 3 ) is used. The mass of 70 Br is not known experimentally but as it in the 
middle of a region of nuclei with measured masses close to the extrapolations of Audi and Wapstra, we have adopted 
their extrapolated value for 70 Br. The mass of 78 Y is not known and not included in the fit. The binding energies 
of the lowest T = 1 states in even-even TV = Z nuclei are derived from those of the isobaric analogues (also taken 
from Ref. [^lj) after an appropriate Coulomb correction. The evolution of the splitting between (J 7r ,T) = (0 + ,l) 
and (J n ,T) = (l + ,0) states in odd-odd nuclei is of particular interest as regards the question of T = 0 and T = 1 
pairing and is currently the object of several experimental studies. The (0 + , 1) state in 58 Cu lies 0.202 MeV above 
the (l+,0) ground state (25). This order is reversed in 62 Ga where the (l + ,0) state is 0.571 MeV above the (0 + ,l) 
ground state f26l| . The E(Tr~, 1) — IT(l + ,0) splitting then continues to rise to 0.837 MeV in 66 As [p7[ . A very recent 
experiment on™Br ^8| has not observed a (l + ,0) level; the lowest observed T = 0 level (with J 77 = 3 + ) is at an 
excitation energy of 1.337 MeV. Similarly, the lowest T = 0 state in 74 Rb measured by Rudolph et al. [GO) at an 
excitation energy of 1.006 MeV has J = 3 and the energy of the 1 + state is unknown. With these data as input, the 
parameters in (@) can be adjusted through a fit procedure that minimizes the rms deviation in the binding energies of 
two states per nucleus (if known). The resulting parameters are shown in the line labeled t56 Ni to 78 Y’ of Table || and 
lead to an rms deviation of 0.396 MeV. In Fig. the differences in energy between the T = 1 and T = 0 states are 
compared to the observed ones. One notes the good agreement that is obtained which gives confidence in the energy 
splittings of 0.847, 1.037, and 1.214 MeV predicted in 70 Br, 74 Rb, and ,8 Y, respectively. As already mentioned, the 
energy difference FT(0 + ,1) — iT(l + ,0) is not known experimentally in these isotopes. In the former two, '°Br and 
74 Rb, the energy difference with the lowest (known) T = 0 state is shown in Fig. |2j. To emphasize the point that these 
energy splittings result from a calculation of total binding energies, the odd-odd results are represented in a different 
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way in Fig. |^. Note that this plot implies a comparison of absolute binding energies: for representation purposes the 
meausured binding energy of the ground state of a particular nucleus is drawn at zero and the other levels of that 
nucleus are given relative to that ground-state energy. 

For the second-half of the 28-50 shell the situation is more complicated since there are no data available. The 
core is 100 Sn with a ground-state mass measured in Ref. [^0|. Since so little is known experimentally, we use the 
extrapolations from Audi and Wapstra ( 2 I]] for the masses of even-even and odd-odd nuclei, complemented with the 
results for ,8 Y from the fit to the first half of the 28-50 shell. The resulting parameters are shown in the line labeled 
t78 Y to 100 Sn’ of Table |. The predictions for the splitting between T = 1 and T = 0 states for the entire 28-50 shell 
are shown in Fig. [|. One notes a satisfactory agreement with the data, when available. The use of extrapolated data, 
however, should weaken the confidence in the predictions for the E{ 0 + , 1) — E(l + , 0) splitting in odd-odd nuclei. 


IV. CONCLUSIONS 

A simple mass formula based on IBM-4 has been proposed to estimate the binding energies of the lowest T = 0 
and T = 1 states of self-conjugate nuclei. It has linear and quadratic terms in the boson number that account for 
the smooth variation of mass with particle number, supplemented with three contributions that have a clear physical 
meaning: an SU(4), a spin-orbit and a T 2 term. It can be considered as a local mass formula that gives predictions 
of a specific interest to current experiments at the N = Z line. As an application we considered nuclei from 56 Ni to 
78 Y where predictions could be made for some of the heavier isotopes currently under study. Also the second half of 
the 28-50 shell was considered although there predictions are more questionable due to the lack of reliable data. 

The advantage with respect to previous IBM-4 work |l7j is that the Hamiltonian used is much simpler and that 
only the L = 0 channel is considered. The numerical diagonalization then becomes trivial and the calculations can be 
performed, without much effort, for arbitrary numbers of bosons. This is much harder to achieve with the full version 
of IBM-4. On the down side it should be noted that this approach is restricted to N = Z odd-odd nuclei since other 
odd-odd nuclei have a dominant U(6) representation which is non-symmetric and which cannot be constructed from 
s bosons only. Also, deformation effects which are present with s and d bosons and which can be included through 
orbital operators [associated with algebras represented by dots in (jl])], are outside the scope of the simple approach 
presented here. 
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TABLE I. Parameters (in MeV) for the 8-20 and 28-50 shells. 


Shell 

a 

P 

7 

£ 

V 

lb O to 28 Si 

16.060 

0.477 

0.190 

-6.146 

-3.009 

30 P to 40 Ca 

-24.538 

0.110 

0.065 

-3.735 

-1.846 

56 Ni to 78 Y 

22.815 

0.118 

-0.068 

-1.953 

-0.450 

78 Y to 100 Sn 

-28.464 

0.118 

-0.188 

-1.045 

-0.512 


TABLE II. Binding energies (in MeV) of IV = 
the parameters given in Table jlj. 

Z nuclei in the first half of the sd shell. Calculated values 

are obtained with 

Nucleus 

T 

-fi'Expt 

ElBM4 

A E 

rap 

0 

151.662 

152.573 

-0.912 

ra p 

1 

150.620 

152.701 

-2.081 

20 Ne 

0 

178.307 

178.887 

-0.580 

20 Ne 

1 

168.033 

167.755 

0.278 

22 Na 

0 

195.476 

195.332 

0.143 

22 Na 

1 

194.819 

194.722 

0.097 

24 Mg 

0 

223.545 

222.918 

0.628 

24 Mg 

1 

214.029 

212.543 

1.486 

26 A1 

0 

241.423 

242.181 

-0.758 

26 A1 

1 

241.195 

240.774 

0.421 

28 Si 

0 

270.581 

271.029 

-0.448 

28 Si 

1 

261.265 

261.465 

-0.200 


TABLE III. Binding energies (in MeV) of N = 
with the parameters given in Table jlj. 

Z nuclei in the second half of the sd shell. Calculated values 

are obtained 

Nucleus 

T 

-S'Expt 

-ElBM4 

A E 

3Up 

0 

289.433 

289.456 

-0.024 

30p 

1 

288.756 

288.968 

-0.213 

32 S 

0 

315.655 

315.300 

0.350 

32 s 

1 

308.653 

308.507 

0.146 

34 Cl 

1 

334.744 

334.723 

0.021 

34 Cl 

0 

334.598 

334.938 

-0.340 

36 Ar 

0 

361.450 

361.513 

-0.063 

36 Ar 

1 

354.839 

354.456 

0.383 

38 k 

0 

381.186 

381.351 

-0.165 

38 k 

1 

381.056 

381.393 

-0.337 
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FIG. 1. Calculated energy differences E(T = 1) — E(T = 0) in N = Z sd-shell nuclei for the parameters given in Table | 
compared with the experimental differences and those of Macchiavelli ef al. [E2|. 
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FIG. 2. Calculated energy differences E(T = 1) — E(T = 0) in N = Z nuclei between 58 Cu to 78 Y for the parameters given 
in Table |J, compared with the experimental differences and those of Macchiavelli et al. ^2j. 
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